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We evaluate the nucleon axial form factor, GA(g^), and induced pseudoscalar form factor, Gp(q^), 
as well as the pion- nucleon form factor, G,r]vjv(g^), in lattice QCD. We also evaluate the correspond- 
ing nucleon to A transition form factors, C^{q^) and C^{q^), and the pion-nucleon-A form factor 
GnNAiq'^). The nucleon form factors are evaluated in the quenched theory and with two degenerate 
flavors of dynamical Wilson fermions. The nucleon to A form factors, besides Wilson fermions, 
are evaluated using domain wall valence fermions with staggered sea quark configurations for pion 
masses as low as about 350 MeV. Using these form factors, together with an evaluation of the 
renormalized quark mass, we investigate the validity of the diagonal and non-diagonal Goldberger- 
Treiman relations. The ratios Gttna {q^)/G^NN{q^) and 2Ci{q^) /GA{q^) are constant as a function 
of the momentum transfer squared and show almost no dependence on the quark mass. We confirm 
equality of these two ratios consistent with the Goldberger-Treiman relations extracting a mean 
value of 1.61(2). 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 



I. INTRODUCTION 

Form factors measured in electromagnetic and weak 
processes are fundamental probes of hadron structure. 
Despite the long history of experimental [IJ and theoret- 
ical studies[2] on nucleon electromagnetic form factors, 
new measurements of these quantities continue to reveal 
interesting features. The discrepancy between the ra- 
tio of the electric to magnetic nucleon form factors ex- 
tracted via Rosenbluth separation and from recent po- 
larization measurements is a well known example. The 
transition form factors in jN A have recently been 
measured jS), |4| to high accuracy, paving the way for theo- 
retical studies using chiral effective theories [a| and lat- 
tice QCD 7, 8, 9] . Compared to the electromagnetic form 
factors, the nucleon (N) and nucleon to A form factors 
connected to the axial- vector current are more difficult to 
measure and therefore less accurately known. An excep- 
tion is the nucleon axial charge gA = C?a(0), which can be 
determined precisely from /?— decay. Its g^-dependence 
has been studied from neutrino scattering [lO| or pion 
electroproduction [ll|, On the other hand the nu- 

cleon induced pseudoscalar form factor, Gp{q^), is less 
well known. Muon capture and radiative muon capture 
are the main experimental sources of information p^ . 
Both GA{q^) and G„{q'^) have been discussed within chi- 
ral effective theories [H, [l^. The electroweak N to 
A transition form factors are even less studied. Using 
Adler's parametrization [l^ the N to A matrix element 
of the axial vector current can be written in terms of four 
form factors, two of which are suppressed fl7]. The two 
dominant transition form factors, C^{q'^) and C|^(g^), 
are analogous to GA{q^) and Gpiq"^), respectively. Neu- 
trino interactions in hydrogen and deuterium were stud- 



ied in an effort to extract information on these form 
factors. Experiments using electroproduction of the A 
resonance are under way [19| to measure the parity vio- 
lating asymmetry in N to A, connected to leading order 
to the form factor G^{q^). Theoretical input on these 
form factors is therefore very timely and important. 

State-of-the-art lattice QCD calculations can yield 
model independent results on these axial form factors, 
thereby providing direct comparison with experiment. 
Lattice studies reflect the experimental situation regard- 
ing our knowledge of these form factors. There have been 
several recent studies on the electromagnetic nucleon [l^, 
[2l|,[2l| and N to A form factors fl, i, H S, Pi . There have 
also been several lattice evaluations of qa [23, 24, 25|], but 
only very recently lattice studies have began probing the 
g^-dependence of the nucleon axial form factors [la, \2T\ 
and the N to A transition form factors [1^. A notable 
exception is an early lattice study on the nucleon axial 
form factors carried out in th e q uenched approximation 
for rather heavy pion masses [29|. 

In this work we calculate the nucleon axial form fac- 
tors using Wilson fermions in the quenched theory and 
with two degenerate flavors of Wilson fermions [lo, |3l| . 
The lowest pion mass in the case of dynamical Wilson 
fermions that we use is about 380 MeV. We also evaluate 
the pion - nucleon {ttNN) form factor G'^ArAr(g^). For 
the extraction of this form factor we need the renormal- 
ized quark mass, which we calculate via the axial Ward- 
Takahashi identity ( AWI) . In addition, we present results 
on the dominant axial N to A transition from factors 
C^(g2) C^{q^). The pion-nucleon-A (ttTVA) form 
factor, G-^NAiq"^), is also computed in an analogous way 
to the evaluation of G.^NN{q^)- Like in the case of the 
nucleon axial form factors, the starting point is an evalu- 
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ation in the quenched theory using the standard Wilson 
action. A quenched calculation allows an efficient check 
of our lattice techniques by enabling a computation of the 
relevant quantities on a large lattice minimizing finite 
volume effects and obtaining accurate results at small 
momentum transfers reaching pion mass, to^, down to 
about 410 MeV. In the case of the N to A transition, 
the light quark regime is studied in two ways: Besides 
using configurations with two degenerate flavors of dy- 
namical Wilson fermions we use a hybrid combination of 
domain wall valence quarks, which have chiral symmetry 
on the lattice, and MILC configurations generated with 
three flavors of staggered sea quarks using the Asqtad 
improved action [32l |. The effectiveness of this hybrid 
combination has recently been demonstrated in the suc- 
cessful precision calculation of the nucleon axial charge, 
gA as well as in our first evaluation of the N to A ax- 
ial transition form factors ■ Since Wilson fermions 
have discretization errors in the lattice spacing, a, of 0(a) 
and break chiral symmetry whereas the hybrid action 
has discretization errors oi 0{a?) and chirally symmetric 
valence fermions, agreement between calculations using 
these two lattice actions provides a non-trivial check of 
consistency of the lattice results. In this work we obtain 
results on the dominant axial N to A form factors 
and with improved statistics as compared to their 
evaluation in Refs [H, . 

The evaluation of the axial form factors as well as 
the ttNN and ttA^A form factors allows us to check the 
Goldberger-Treiman relations. It is advantageous to cal- 
culate ratios of the nucleon elastic and transition form 
factors, since ratios have weak quark mass dependence 
and we expect them to be less sensitive to other lattice 
artifacts. In particular, it is useful to consider ratios 
for which the renormahzed quark mass cancels since this 
eliminates one source of systematic error. 

The paper is organized as follows: In Section II, we 
give the definition of the matrix elements that we con- 
sider in terms of the form factors on the hadronic level. 
In Section III, we give the lattice matrix elements on the 
quark level and in Section IV, we discuss how we extract 
the form factors from lattice measurements. In Section 
V, we present our results. Finally in Section VI, we sum- 
marize and conclude. 



pseudoscalar currents defined by 



p-{x)^i,{x)^J-^^{x) 
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where are the three Pauli-matrices acting in flavor 
space and if) the isospin doublet quark field. 



A. Axial form factors 

The matrix element of the weak axial vector current 
between nucleon states can be written in the form 
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where we specifically consider the axial isovector current 
A'^. The form factors depend only on the momentum 



transfer squared, = {p' — p^){p 



P^^ 



As defined 



above, the form factors GA{q^) and Gp{q^) are dimen- 
sionless. As already mentioned in the Introduction, there 
exist several recent lattice studies on the nucleon axial 
charge gA [H, [H, [2^ , whereas only very recently there 
are lattice studies to investigate the dependence of 
GA{q^) or Gp{q^), apart from an early calculation in the 
quenched approximation [2^. 

The invariant N to A weak matrix element is expressed 
in terms of four transition form factors (igI . 34 | as 



<A{p',s')\At,\N{p,s) > 

-<A(J2 
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u{p,s) (3) 



where, as in the nucleon case, we consider the physically 
relevant axial isovector current A^^{x). 



B. Pseudoscalar matrix elements 



II. DEFINITION OF MATRIX ELEMENTS 

To extract the form factors, we need to evaluate 
hadronic matrix elements of the form < h'\Ofi\h >, where 
h and h' are the initial and final hadron states and 
a current that couples to a quark. In all that follows, 
we assume isospin symmetry and take the mass of the u 
and d quark to be equal. We consider nucleon-nucleon 
and nucleon-A matrix elements of the axial vector and 



Spontaneous symmetry breaking couples pions to the 
broken axial charges and currents. The relation 

{0\A;{Q)\7t''{p)) = tUp.S'^' (4) 

can be used to extract the pion decay constant, /tt, on 
the lattice by evaluating two point functions. With our 
conventions /tt = 92 MeV. Taking the divergence of the 
axial vector current we obtain the operator relation 

a; = Umln'^ (5) 
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known as the partially conserved axial vector current 
(PCAC) hypothesis. On the QCD level we have the axial 
Ward-Takahashi identity 

(6) 



where all quantities appearing in Eq. ([6]) are renormalized 
quantities with niq being the renormalized quark mass. 
Comparing Eqs. ([5]) and ^ we can write the pion field 
7r° in terms of the pseudoscalar density as 



(7) 



The renormalized quark mass can be evaluated by taking 
the matrix element of Eq. Q between a zero momentum 
pion state and the vacuum to obtain 



< 0|Ag|7r°(0) > 
2 < 0|P''|7r'^(0) > 



(8) 



Taking the matrix clement of the pseudoscalar density 
between nucleon states we can define the ttNN form fac- 
tor via 
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2mg < N{p')\P^\N{p) > = 



ip')ij5u{p) ■ (9) 



Similarly the nuclcon-A matrix element of the pseu- 
doscalar density yields the TriVA form factor: 



27nq < A{p')\P''\N{p) >= 



EA{p')ENip) 
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2 2 ^'^(P'o • (10) 

Eqs. ^ and pH]) define the form factors G-^NNiq'^) 
and GTrNAiq"^) that we study in this work. The ttNN 
and TriVA strong coupling constants are then given by 
g^NN = Gt,nn{0) and g^NA = G^jva(O). PCAC re- 
lates the axial form factors Ga and Gp with Gt^nn and 
equivalently and Cq with Gt^na- Using the PCAC 
hypothesis together with Eq. ^ we obtain the diagonal 
Goldberger-Treiman relation (GTR) 



G,(,^) + -^G,(,^)^^^^^^#^ .(11) 

4m^ 



2mN 
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Similarly using Eq. (fTO)) we obtain the non-diagonal GTR 
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Assuming pion pole dominance we can relate the form 
factors Gp to Gt^nn and to G-^na via: 

2G^Ar7v((7^)/7r 



1 

2mN 
1 

mjv 



C^{q') 



ml — q^ 

IGjrNAijf)]^ 

2 ml — q^ 



(13) 



Substituting in Eqs. pTjl and we obtain the simpli- 
fied Goldberger-Treiman relations 



GjrNNiq^) /tt 

G^ArA('Z^) /tt 



mNGA{q^) 
2mNG^{q^) 



(14) 



III. LATTICE EVALUATION OF 
CORRELATION FUNCTIONS 

To evaluate the axial nucleon form factors Ga and Gp, 
we use the techniques developed in our study of the nu- 
cleon isovector electromagnetic form factors [l^]. Since 
only the isovector axial vector is of relevance here, only 
the connected diagram shown in Fig. [T] is needed. To 
extract the matrix element of the axial isovector current 
between nucleon states defined in Eq. ([2]) we need to cal- 
culate the three-point function 

(G^^>(t2,ti;p',p;r)) 
= ^ exp(-ip ' ■ X2) exp(+i(p ' - p) ■ xi) r''" 

(J7|T[x"(x2,t2)4(xi,ti)x''(0,0)] 1^7), (15) 

using the local quark bilinear axial current A^^(x) of 
Eq. ID). An interpolating field with the quantum num- 




FIG. 1: Connected three-point function between final and 
initial hadron states /i'(p') and h{p). 

bers of the nucleon routinely used in lattice studies is 

X{x)=e''^'' [u^''{x)C-izd\x)\ u'^ix) . (16) 

We create an initial state (source) by acting with x(0) 
on the vacuum. Evolution in Euclidean time with the 
QCD Hamiltonian produces, in the large time limit, the 
nucleon state. We take ti, the time from the source at 
which the axial-vector current couples to a quark, to be 
large enough so that the nucleon is the dominant state. 
We then take the overlap with a nucleon state that is an- 
nihilated at a later time ^2 — by the same interpolating 
field x{x) (sink). Again we take ^2 ^ large enough so 
that the nucleon is the dominant state. This process is 
schematically shown in Fig. [TJ In addition, we calculate 
the nucleon two-point function. 



(G'^^(t,p;r)>=^, 



{n\Tx"(^,t)x''{o,o) |n>(i7) 
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where the projection matrices for the Dirac indices are 
given by 



sidered is 



cr,; 




/ 




(18) 



We then construct a ratio such that, in the large Eu- 
chdean time hmit all exponential dependences on times 
and unknown initial state-nucleon overlap constants < 
A^|x|0 > cancel. One can find more than one ratio that 
accomplishes this. We require, in addition, that we use 
two point functions that involve the shortest time evolu- 
tion. Such a ratio is given by 



'(t2,ii;p',p;r)> 



(G^^(t2-ti,p;r4)) (G^^(ti,p';r4)) (G^^(t2,p';r4)> T^/' 



{G^^(t2-ti,p';r4)) (G^^(ti,p;r4)) (G^^(t2,p;r4)> 
=> n (p ,p ;r;/i) , 



(19) 



which, in the large Euclidean time, with the nucleon 
the dominant state, produces a constant plateau region 
in ^1. Throughout this work we use kinematics where 
the final hadron state is produced at rest and therefore 
the momentum transfer q = p' — p = —p. We take 
—q^ = > with being the Euclidean momentum 
transfer squared. The value of R'^{t2, ii; p ', p ; F; fi) in 
the plateau region, n'^(p', p ; F; /z), is directly connected 
to the nucleon form factors through the relation 



n^(0,-q;rfc;A*) 



C 



{{En + mN)5k,^ + qk5^A)GAiQ^) 



2m 



N 



Gp(Q2) 



(20) 



for k = 1,2,3, while 0-^(0, -q \Ti;p) = 0. The nu- 



2mi 



cleon energy En = ^mj^ + q2 and C = ^ E«(E«+mM) ' 
a factor related to the normalization of the lattice states. 
Since our goal is to evaluate the form factors as a func- 
tion of , we calculate the three point functions with se- 
quential inversions through the sink. This requires fixing 
the source-sink time separation t2 as well as the initial 
and final hadron states but allows the insertion of any 
operator with arbitrary momentum and any time slice 
ti. In fact, the usefulness of this technique is evident in 
this work: Since the same matrix elements were calcu- 
lated for the electromagnetic current [1, no new 
sequential inversions are required for the axial current or 
pseudoscalar density operators. 

As in the electromagnetic case, it is advantageous to 
use a linear combination of nucleon interpolating fields 
to construct optimal sources and sinks. Since for axial 
operators a non-zero contribution can be obtained only if 
r ^ r4 in the three-point function, the most symmetric 
linear combination of matrix elements that can be con- 



5-"(q; j)-En-'(o.-q;rfc;/^ = j) = 



k=l 



c 



4m 



N 



{En + mjv) ((5ij + (Jsj + h^j) Ga{Q^) 
-{qi+q2 + q3)^Gp{Q^) 



(21) 



where j — 1,2,3 labels the spatial current direction. In 
order to obtain the three matrix elements corresponding 
to three choices of Ffc in the above sum one would require 
three sequential inversions. However, choosing an appro- 
priate linear combination of nucleon interpolating fields, 
this sum is automatically built in and with one sequential 
inversion we can obtain S''^(q; j) for all current directions 
j. We call such a linear combination optimal sink because 
it allows us to take into account in our determination of 
the form factors the largest set of momentum vectors 
contributing to the same value. Since the sequential 
propagators corresponding to this sink have been com- 
puted for the isovector electromagnetic form factors [2^ 
they can be used directly here. Therefore the computa- 
tional cost for the evaluation of the three-point function 
for all intermediate times ti, current indices /i and a large 
set of lattice momenta vectors q is very small. 

Similarly, to evaluate the ttNN form factor G-^nn we 
construct the ratio i?^, which is the same as the ratio 
given in Eq. but instead of the three-point func- 
tion (G'^^'^(i2,ii;p',p;r)) defined in Eq. ^ with 
the axial current A'^^, we use the three-point function 

(G^-f''^(t2,ii;p',P;r)), obtained by replacing in 
Eq. 1131) by the pseudoscalar density . The large Eu- 
clidean time behavior of is independent of ti leading 
to the plateau value denoted by n'^(0,— q ;r;75). The 
value of R^ in the plateau region, 11-^, is related to the 
ttNN form factor via 

n^(0,-q;rfe;75) = 

C 7 — 2 'I G-knn{Q ) (22) 

2m N 2mq{m^ + Q^) 

for A: = 1,2,3 while 11^(0, —q ;r4;75) = 0. Summation 
over the polarized matrix elements now leads to 



^''(q; 75) = ^n^(o,-q;rfe;75) = 

k=l 

^ 2mN 2m,(m^ + Q2)G.A^A'W) '(23) 



from which G-^nn can be extracted if we know f.^ and 
mq. 

The determination of the N to A axial transition form 
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factors requires the evaluation of the three-point function 

(G„ " (t2,ti;p',p;r)) 
= ^ exp(-ip ' ■ X2) exp(+i(p ' - p) ■ xi) r'''^ 

{n|r[x"(x2,i2)4(xi,ti)x''(0,0)] |n> ,(24) 

where to create an initial state with the A+ quantum 
numbers we use the standard Rarita-Schwinger interpo- 
lating field 

+ [u^''{x)Cj„u''{x)]d'{x)Y (25) 

Besides the nucleon two-point function we also need the 
A two-point function given by 

= ^e-p'-^r"" (nirx"(x,t)x?(o,o) Ifi) . (26) 

X 

The corresponding ratio which, in the large Euclidean 
time limit becomes ti-independent, is given by 



R^{t2,ti;p',p ;r;/i) 



{G^^(t2,p';r4)) 



(G,^^(t2,p';r4)) (G'^^(t2-ti,p;r4)> (G^'^(fi,p';r4)) 



(G'^^(t2,p;r4)) (G^^(t2-ti,p';r4)> {G^^(fi,p;r4)) 
=> (p ,p ;r;/i) . 



1/2 



(27) 



As in the nucleon case, we take the final A state to be 
produced at rest and therefore q = p' — p = — P- The 
value of in the plateau region, U^, for the case T = 
which is the relevant one for this work, is related to 
the form factors via 



n^(0,-q;r4;j) 



iB 



{En - 2mA + ^'^a') 



2{EN + mN) 





niN 



(28) 



for the spatial components of the current, /i — j, whereas 
for the temporal current component, /i = 4, we have 



n^(0,-q;r4;4) = V 



mA ^ -Ejv - TOA 



(29) 

for k — 1,2,3, while 114(0, —q; F; /i) — 0. In this case 
we have a larger freedom in choosing appropriate linear 
combinations for the optimal sink due to the additional 
vector index of the A. Using this freedom we construct 
A sinks so that the maximum allowed number of lattice 
vectors contribute in the evaluation of the form factors 
at a given value of Q^. These turn out to be the same 



as the ones used in our studvof the N to A electromag- 
netic transition form factors [8| . Therefore the sequential 
inversions already performed for the evaluation of the N 
to A electromagnetic transition form factors d, Q can 
be used for the computation of the axial transition form 
factors. We give below the expressions that we obtain in 
the large Euclidean time limit, using the optimal linear 
combinations of A interpolating fields: 



Si{ci;j) = ^n^(0,-q;r4;j) =jB 



J {En - 2mA + mjv) + ( ^P^' 



ruA 

TUN 



{En - mA)Ct + ruNCi - ^^p' { Y^p 

ruN ' ^-^ 



En + TUN 
3 ^ 



,(30) 



Si^(q;4) = ^n^(0,-q;r4;4) = i3 

o- = l 



fc=l 

3 



, rriA „A , En - „a 

O3 H O4 H Oe 

mjv mjv 



(31) 



S2^(q;j) 



n^(0,-q;rfc;i) = 



.34 



- p') + SjAp"" - p') + SjMp' - p')) 



S^{ci;j) = n^(0,-q;r3;M)-^ 



nf(0,-q;ri;^i) 



-n^(0,-q;r2;M) 



= iA 



fx 2 3 r 1 3\ 

(Oj.ip p - dj,2P p j 63 



(32) 



,(33) 



where j = 1,2,3 and 



A 



B 



{Em + niN) 



B = 



2 ^{En + niN) /En 



3m 



N 



(34) 



As can be seen, S'^(q; j) and S^{q;j) isolate the sup- 
pressed form factor for different combinations of lat- 
tice momentum vectors as compared to that also in- 
volves C^. Since here we are only interested in the dom- 
inant form factors and Cq, we use only Si{q;fi). 
We denote the ratio constructed analogously to R^, but 
with the optimal sink, by Rna- T^^^ linear combi- 
nation (q; /i) turns out to be the suitable one also 
for the calculation of the from factor G^ina, defined 
in Eq. pO]) . Again replacing the three-point function 

{Go- " (t2, ii; p ', p; r)) with the corresponding pseu- 
doscalar three-point function {G^^ ^ {t2,ti;p ' ,p;T)) in 
Eq. ([7f|) . we obtain the ratio R^ , which at large Eu- 
clidean times becomes time-independent. Fitting in the 
plateau region yields 11^(0,— q ;r;75), related to the 



6 



TriVA form factor via the relation 



pion state with the vacuum we form the ratio 



n^(0,-q;r4;75) 



2 /En + niN Qo 



En Gtun '2mq{ml + Q^) 
The optimal combination gives 



^^A(q; 75) = 5]n^(o,-q;r4;75 

(7 = 1 

qi+q2 + qs firinl 



(35) 

2 I r)2A '^^na{Q ) 



2 / En + TUN 
3' 



En 

G^NAiQ^). (36) 



from which Gj^na can be determined. 

In order to evaluate G^at at and G^tva , we need to know 
niq . The renormalized quark mass can be defined via the 
AWI given in Eq. ([8]). On the lattice the AWI has correc- 
tions, which in the case of Wilson fermions are of order a. 
For domain wall fermions (DWF), the divergence of the 
four-dimensional vector axial current has an additional 
term that goes to zero as the fifth dimension goes to infin- 
ity [ssj . For non-singlet matrix elements at low energies, 
this additional term shifts the quark mass by an additive 
constant known as the residual quark mass [s^. Such a 
simple shift is valid up to order a^. Provided that the 
fifth dimension is large enough, the residual mass maybe 
considered a small correction. The size of the fifth di- 
mension was adjusted by requiring the residual mass to 
be small compared to the pion mass. The criterion used 
to fix the fifth dimension to 16 is that the residual mass 
is smaller than 10% of the quark mass [37j. We adopt 
the same criterion in this work and take the fifth di- 
mension to be 16. Therefore we assume that corrections 
due to the residual mass are small and we calculate the 
renormalized quark mass using the AWI taking matrix 
elements between a pion zero momentum state and the 
vacuum. The initial state with the pion quantum num- 
bers is created using the axial-vector current as an 
interpolating field. The quantities with the tilde denote 
operators that are built from smeared quark fields ob- 
tained from point quark fields ip(x) as described in the 
next subsection. The pion-vacuum matrix element of the 
axial-vector current is given by the two-point function 

CUt) = E ^ [Al{^,t)Al{0,0)) \n ) (37) 

X 

and the pion-vacuum matrix element of the pseudoscalar 
density is given by 



The subscript L denotes that the axial- vector current and 
pseudoscalar density are constructed using local quark 
fields unlike the interpolating fields A^l and that use 
smeared quark fields . To cancel the overlaps of our initial 



2 ZpC£s{t)\ Cisit) 



(39) 



using, in addition to local-smeared (LS) two-point func- 
tions, the smeared-smeared two-point functions Cgg and 
Cgg. We look for a plateau in the large Euclidean time 
behavior of the effective mass moff(t), which determines 
ruq. The factors Za and Zp are the renormalization con- 
stants for the local axial-vector and pseudoscalar cur- 
rents, respectively. We note that Zp is only needed for 
the determination of the renormalized quark mass. This 
dependence cancels in all physical quantities presented in 
this work, which are therefore independent of the value 
we use for Zp. The evaluation of mg together with the de- 
termination of f-jr allows us to evaluate G-^nn and G^ata • 
The pion decay constant is determined from the large 
Euclidean time behavior of the ratio 



ffit) 



it) 



(40) 



For large the above quantity becomes t independent 
and the plateau value gives Z^. 

For the evaluation of all N to A matrix elements in 
the case of Wilson fermions, we use the sequential prop- 
agators already computed in our study of the electro- 
magnetic transition form factors . However in the 
hybrid scheme, additional propagators are calculated to 
improve the statistical errors beyond those of our previ- 
ous work [9', '28'| and to check for finite volume effects. 
We summarize in Tables III and ITTl the details of the calcu- 
lation. All the hadron masses given in Table [Til are com- 
puted using domain wall valence quarks and MILC con- 
figurations for the sea quarks. The value of the valence 
domain wall quark mass, , was determined by tun- 

ing the pion mass calculated with domain wall fermions 
to be the same as the lowest mass pion in the staggered 
formulation [ST*]. 

For Wilson fermions, we convert dimcnsionless lattice 
quantities to physical units by setting the lattice spacing 
using the nucleon mass at the physical limit. The value 
of a extracted from the nucleon mass is given in Table U 
and it is consistent with the value extracted using the 
Sommer scale rg. The dynamical Wilson configurations 
at K = 0.1575 and 0.1580 were generated by the TxL col- 
laboration [sSl and at K = 0.15825 by the DESY-Zeuthen 
group [Slj. For the hybrid calculation we use the scale 
extracted from heavy meson spectroscopy as determined 
by the MILC collaboration j33|- As can be seen in Ta- 
ble ini in the hybrid approach we consider lattices with 
temporal extent 32 and 64. Temporal extent 32 is ob- 
tained by using Dirichlet boundary conditions (b.c.) in 
the temporal direction to cut into half the original MILC 
lattices when we calculate the domain wall quark prop- 
agator. This was the procedure adopted in our previous 
evaluation of N to A axial form factors Q due to the 
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Wilson fermions 



No. of confs K m,r (GeV) m]v (GeV) tua (GeV) 

Quenched 32^ x 64, /3 = 6.0, a'^ = 2.14(6) GeV 



200 


0.1554 0.563(4) 


1.267(11) 


1.470(15) 


200 


0.1558 0.490(4) 


1.190(13) 


1.425(16) 


200 


0.1562 0.411(4) 


1.109(13) 


1.382(19) 




=0.1571 0. 


0.938(9) 




Unquenched[30] 24^ x 40,/3 = 5.6, 


= 2.56(10) GeV 


185 


0.1575 0.691(8) 


1.485(18) 


1.687(15) 


157 


0.1580 0.509(8) 


1.280(26) 


1.559(19) 



Unquenched[31] 24^ x 32,/3 = 5.6, = 2.56(10) GeV 



200 0.15825 0.384(8) 1.083(18) 1.395(18) 
= 0.1585 0. 0.938(33) 



TABLE I: Parameters for the calculations using Wilson fermions 







Hybrid action a 


-1 = 1.58 GeV [38] 






No. of confs 


Volume 


(amu,d) 


/ \ sea 

[anis ) 


/ \D\V 

(amq) 


(GeV) 


TUN (GeV) 


ruA (GeV) 


150 


20^ x 32 


0.03 


0.05 


0.0478 


0.606(2) 


1.392(9) 


1.670(22) 


150 


20^ X 32 


0.02 


0.05 


0.0313 


0.502(4) 


1.255(19) 


1.567(25) 


118 


28^ X 32 


0.01 


0.05 


0.0138 


0.364(1) 


1.196(25) 


1.561(41) 


200 


20^ X 64 


0.03 


0.05 


0.0478 


0.594(1) 


1.416(20) 


1.683(22) 


198 


20^ X 64 


0.02 


0.05 


0.0313 


0.498(3) 


1.261(17) 


1.589(35) 


100 


20^ X 64 


0.01 


0.05 


0.0138 


0.362(5) 


1.139(25) 


1.488(71) 


150 


28^ X 64 


0.01 


0.05 


0.0138 


0.357(2) 


1.210(24) 


1.514(41) 



TABLE IL Parameters for the calculationz using the hybrid action 



limited computer resomxes. In this work we present, in 
addition, results using the full temporal extent of the 
MILC lattices with antiperiodic b.c. in the temporal di- 

I 



IV. EXTRACTION OF OBSERVABLES 

A. Ground state dominance and noise reduction 

As we already pointed out, in order to extract physical 
matrix elements, we must first evolve in Euclidean time 
to create the hadronic state of interest. In this work, 
the hadronic states of interest are the pion, the nucleon 
and the A states. To create the initial states with the 
pion quantum numbers, we use the temporal component 
of the axial- vector current, and for the nucleon and A, 
the interpolating fields given in Eqs. and re- 
spectively. The length of the time evolution required to 
obtain the true pion, nucleon and A eigenstates depends 
on our choice of the initial state. It is well known that 
if one constructs a hadron initial state using smeared 
quark fields instead of localized ones, the convergence to 



rection consistent with what is used in the simulation 
of the configurations. Antiperiodic b.c. in the temporal 
direction are also used in the case of Wilson fermions. 



I 

the hadron cigcnstate is very much improved. Therefore, 
in this work, we always smear the quark fields in a gauge 
invariant way using the so called Wuppertal or Gaussian 
smearing |39j . In this scheme the smeared quark field, 
tp{x,t), is obtained from the localized field, ip{z,t), via 

^(x,i) =^F(x,z;[/(t))^(z,i) . (41) 

z 

The gauge invariant smearing function is constructed 
from the hopping matrix H: 

Fi^, z; U{t)) = (1 + ai7)"(x, z; U{t)), (42) 

where 

i7(x, Z; Uit)) = J2 t)S^,y-^ + U}{^- i, t)S^^y+^^ . 

i=l ^ ^ 

(43) 
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The parameters for the Wuppertal smearing are deter- 
mined by requiring that the nucleon state dominates the 
two-point correlator for the shortest time evolution. We 
find that a = 4 and n = 50 are optimal parameters. 
Although smearing improves ground state dominance, it 
introduces gauge noise increasing the errors on the ex- 
tracted effective masses in particular when applied both 
to the source and to the sink. An efficient way to re- 
duce the ultraviolet fluctuations is to smooth the gauge 
fields that enter the hopping matrix H{x, z; U ( t) ) . It 
was shown in Ref. that hypercubic smearing ji^l on 
these gauge links reduces gauge noise and tends to also 
improve ground state dominance. In the case of domain 
wall fermions, HYP smearing is in fact applied to all 
the gauge links so as to accelerate the convergence of 
the bi-conjugate gradient method used to evaluate the 
inverse of the fermionic matrix. In the quenched case, 
HYP smearing is not used because self averaging is more 
effective on larger lattices. In the case of dynamical Wil- 
son fermions, the simulations are done on smaller lattices 
causing gauge noise, and HYP smearing needs to be ap- 
plied to the gauge fields that enter the hopping matrix 
i?(x,z;[/(i)). 
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B. Plateaus and Overconstrained Analysis 

In this subsection we describe the analysis of the lat- 
tice measurements that lead to the extraction of physical 
quantities. 

The mass of the lowest hadron state for a given set of 
quantum numbers is the simplest quantity to calculate 
on the lattice, since it requires only the computation of 
two-point functions. In this work we need, besides the 
pion, the nucleon and the A mass, which are straightfor- 
ward to determine, the renormalized quark mass. This 
is evaluated by taking matrix elements of the AWI as 
discussed in Section HI. The effective mass m^^(i) de- 
fined in Eq. ([55)1 becomes time-independent if t is large 
enough so that the pion ground state dominates (plateau 
region). We show in Fig. [5] m^^^{t) as a function of 
time, both in lattice units. We consider all three values 
of the bare quark mass for each of the three types of sim- 
ulations that we use in this work, namely the quenched 
approximation, two dynamical Wilson fermions and the 
hybrid scheme. As can be seen, in all cases, allowing for 
an initial time evolution, the effective mass becomes time 
independent. Fitting to a constant in the plateau region 
determines irig. In the case of DWF, the extracted value 
is expected to be the same as used in the domain 

wall Dirac matrix and given in Table [TTl Any differences 
are attributed to an additive residual mass that provides 
a measure of the chiral symmetry breaking due to the 
finite extent of the fifth dimension. In the case of Wil- 
son fermions it is known that the AWI has corrections of 
0{a). The value of therefore carries systematic er- 
rors, which can be large as we approach the chiral limit, 
since corrections that appear in the right hand side of 



t/a 



FIG. 2: The effective quark mass (i) defined in Eq. ^ 

as a function of time, both in lattice units. The upper graph is 
for the quenched case, the middle graph for dynamical Wilson 
fermions and the lower graph for the hybrid scheme. The 
dashed lines span the range of fitted points and show the 



Eq. ^ will dominate as the first term decreases. The 
results on will be discussed in the next section. 

For the evaluation of form factors we look for time 
independence of the ratios R^, , Rna ^^"^ -^^A 
structed from three-point functions and appropriate com- 
binations of two-point functions. We start with the ratio 
R^ , which in the large Euclidean limit goes to given 
in Eq. (|23p . from which the ttNN form factor is deter- 
mined. We show in Fig. [3] typical examples of the ratio 
R^ divided by C{qi + q2 + q3)/2mN and averaged over 
all momentum directions that lead to the same value. 
Such an averaging is also done in our overconstrained 
analysis described below. The ratio is shown as a func- 
tion of ti in physical units for the four lowest Q^-values 
both in the quenched theory at the intermediate quark 
mass and for dynamical Wilson fermions at the heaviest 
quark mass. When the time separation from the source 
and sink is large enough so that the nucleon state dom- 
inates, this averaged ratio is time independent. When 
this happens the quantity plotted in Fig. [3] corresponds 
to fyrmlGTrNN /[2'mq{Q'^ + 7n^)Zp]. Note that, since to 
obtain the renormalized mass we divide by Zp, the Zp 
factors cancel and therefore we do not need to know Zp. 
As already pointed out, this is true for all physical quan- 
tities that we calculate in this work. The dashed lines 
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2 3 

t, (GeV-0 

FIG. 3: The ratio used to extract G^nn for the four lowest 
values of Q^. The upper graph is for the quenched theory 
at K = 0.1558 (m,r ~ 0.49 GeV) and the lower graph for 
dynamical Wilson fermions at k = 0.1575 (m,r ~ 0.69 GeV). 
The dash lines are fits to the plateaus and span the range of 
fitted points. 



show both the range used for the fit and the value of the 
plateau. In Fig. [3]we show the corresponding average of 
the ratio that, in the large time limit, leads to 
defined in Eq. ([55)1 and to the determination of the ttTVA 
form factor. Here we show results for pions of mass of 
about 500 MeV (intermediate value) in each of the three 
types of simulations considered in this work. As can be 
seen the quality of the plateaus in all cases is good enough 
to allow us to fit to a constant within the plateau range 

with a x^/degrees of freedom (d.o.f) ^ 1. This leads to a 
good determination of the ttNN and ttA^A form factors. 

In Fig.[5]we show the ratio defined in Eq.[Tn]for the 
optimal combination that, for large ti and t2 — ti time 
intervals, leads to Eq. pT|) . We show results for the four 
lowest values of the lattice momentum vector, q, since, in 
this case, the ratio depends on the momentum vector and 
not just its magnitude. We note that this is not what is 
actually fitted, since in our overconstrained analysis we 
consider all lattice vectors q that result in the same 
value. However Fig. Ogives an idea of the quality of the 
plateaus that are used in the overconstrained analysis to 
extract the nucleon axial form factors. Similar plateaus 
are obtained for the ratio R'^j^ needed to extract the N to 
A axial form factors. Again the quality of the data allows 
identification of a plateau region to which to perform a 
fit to a constant to extract the matrix element we are 
interested in. 

The overconstrained analysis uses all the stochastically 
independent lattice measurements that contribute at a 
given when extracting the form factors 4l|. This is 
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FIG. 4: The ratio -Rjva used to extract GnNA for the four 
lowest values of Q^. The upper graph is for the quenched 
theory, the middle graph for dynamical Wilson fermions and 
the lower graph for hybrid scheme for a pion of mass about 
500 MeV (intermediate value). The dash lines are fits to the 
plateaus and span the range of fitted points. 



done by solving the overcomplete set of equations 

Pi^;^l)^Diq■^l)■F{Q^) (44) 

where -P(q; fJ-) are lattice measurements of appropriately 
defined ratios. For concreteness, let us consider the anal- 
ysis for the nucleon axial form factors. In this case 
P(q; /i) is the ratio R^ given in Eq. having statistical 
errors Wk- The vector F contains the form factors: 



(45) 



If N is the number of current directions and momentum 
vectors contributing to a given then 13 is an x 2 
matrix which depends on kinematical factors. We extract 
the form factors by minimizing 



Pk 



Wk 



(46) 



using the singular value decomposition of D. Therefore 
we do not actually fit to the plateaus shown in Fig. [5] for 
each momentum vector but combine all momentum vec- 
tors in the overconstrained analysis. A similar analysis is 
done for the determination of all the other form factors. 
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FIG. 5: The ratio defined in Eq. ^ with the opti- 
mal nucleon sink from which Ga and Gp are extracted 
for the four lowest momentum vectors q = (1, 0, 0)27r/Ls, 
q = (1, 1, 0)2n/Ls, q = (1, 1, l)2n/L, and q = (2, 0, 0)2tt/Ls, 
where La is the spatial size of the lattice. The upper graph 
is for the quenched theory at intermediate pion mass (k = 
0.1558) and the lower for two dynamical Wilson quarks at 
the heaviest mass (k = 0.1575). The dash lines are fits to the 
plateaus and span the range of fitted points. 



FIG. 6: The upper graph shows R , from which Ga and 
Gp are extracted, for momentum vectors q — (1, 0, 0)27r/Ls 
and q — {l,l,l)2n /Ls- The lower graph shows the ratio 
, from which G-nNN is determined, for momentum transfer 
squared q^ — {2-k/Ls)'^ and q^ — 2(27r/Ls)^. Results on 
these quantities are shown as a function of t\ for the quenched 
theory for sink-source separations, 12/ a — 13 (crosses), and 
t2/a = 11 (open triangles) at the smallest quark mass (k = 
0.1562). 



C. Fixing the source-sink time separation 

All the results shown in Figs. [H [4] and [5] are obtained 
keeping the source-sink separation, t2, fixed. In the 
quenched case we take a = 11, for dynamical Wilson 
fermions we take ^2/0 — 12 and for the hybrid scheme we 
take t2/a — 8 so as to keep the physical time separation 
approximately constant at about 5 GeV^^ or 1 fm. In 
order to ensure that this time separation is large enough 
to isolate the nucleon and A states we must increase 
the sink-source time separation and check that the re- 
sults remain unchanged. This check is carried out in the 
quenched theory at the lowest quark mass and in the 
hybrid scheme. In both, we increase the source-sink sep- 
aration by two time slices. 

We choose to do this check for quenched rather than 
dynamical Wilson fermions since the errors are smaller 
and we can therefore identify deviations more easily. We 
choose the smallest mass, since the smaller the mass 
the more severe is the contamination of excited states. 
In Fig. [S] we show the ratio for the optimal nu- 
cleon source of Eq. (PT|) for source-sink time sepa- 
ration t2/a = 11 and 12/0, = 13. Results are shown 
for the lowest momentum vector q = {1,0,0)2tt/ Lg and 
for q ~ (1, 1, l)27r/Ls, where Lg is the spatial extent 
of the lattice. As can be seen the ratios yield consis- 
tent plateaus. In the same figure we also show the ratio 



from which G-kNN is extracted for the two lowest 
q^-values. The only discrepancy arises at the lowest q^ 
value, where the larger source-sink separation prodices 
a larger result. At all higher values of q^, the plateaus 
are however consistent as demonstrated for the second 
lowest value of q^. Given that the plateaus for R^ for 
both time separations are consistent at all values of the 
momentum vectors, the discrepancy seen in the case of 
R^ at the smallest q^ value may have a different origin. 
We note that is proportional io qi + q2 + qz- As q 
extracting G-^nn becomes ill-defined and our statistical 
error in this case underestimates the true error. The ef- 
fect of increasing t2 on the actual form factors can be 
seen in Fig. [71 where we show the nucleon axial form fac- 
tors Ga{Q^) and Gp[Q'^) as well as G-j^NNiQ^) extracted 
for sink-source separations 12/ a = 11 and ^2/0 — 13. As 
can be seen, the results up to ~ 1.5 GeV^ at the two 
time separations are within error bars with the only ex- 
ception the value of G^^nn at the lowest value, which 
differs by about one standard deviation. Differences by 
about one standard deviation in the results for GaIQ^) 
for > 1.5 GeV^ are, most likely, due to taking nu- 
merically the Fourier transform, which for large values of 
Q^, becomes noisy, requiring more statistics. Given this 
level of agreement at the smallest quark mass we con- 
clude that, for the quenched physical time dis- 
tance of about 5 GeV~^ ~ 1 fm is sufficient for ground 
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FIG. 7: The upper graph shows G^{Q^), the middle graph 
Gp(Q^) and the lower graph f-nG-nNN /rnq as a function of 
in the quenched theory for sink-source separations, 1-2/ a = 13 
(crosses), and 12/ a — 11 at the smallest quark mass {k — 
0.1562). 
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FIG. 8: The upper graph shows C^{Q^), the middle graph 
C^(Q^) and the lower graph f-nG-nNi^lvaq as a function of 
in the hybrid scheme for sink-source separations, 12! a ~ 10 
(crosses), and t2/a = 8 at the smallest quark mass, namely 
mi = 0.01. 



state dominance and identification of a consistent plateau 
region with the hindsight that G^nn at the smallest 
maybe underestimated by about one standard deviation. 

We next discuss the adequacy of the sink-source sepa- 
ration in the hybrid approach. Pion cloud contributions 
are expected to become important for dynamical quarks 
as the quark mass decreases and one must allow a large 
enough time separation for the pion cloud to develop. 
Therefore it is important to ensure that the time sepa- 
ration t2 is large enough for dynamical quarks with the 
smallest mass. The results for the larger time separation 
are obtained using Dirichlet b.c. at the first time slice 
and at the midpoint of the temporal direction cutting in 
half the lattice size whereas for the smaller separation 
antiperiodic b.c. are used. We compare in Fig. [5] the 
N to A form factors extracted for 12/ a — 10 to those 
obtained with sink-source time separation 12/ a = 8. As 
can be seen all the results at the two time separations, 
including Gt^na at the lowest value, are within error 
bars. Given that we use the same number of configu- 
rations for the two time separations it is obvious that 
we have a big advantage for using the smaller separation 
since errors are reduced by more that a factor of two. 
Given the level of agreement at the smallest quark mass 
for both quenched and hybrid results, combined with the 
advantage of smaller statistical errors, we conclude that 



it suffices to take i2 ~ 5 GeV""'^. Therefore all the results 
given in the next section are obtained with this time sep- 
aration. Furthermore results in the hybrid scheme are 
obtained using the full temporal extent of the MILC lat- 
tices with antiperiodic b.c. in the temporal direction. 

D. Volume dependence 

Another potential source of a systematic error is the 
spatial size of our lattices. Given that for the quenched 
case we use a lattice of spatial size of about 3 fm we 
expect finite volume effects to be negligible. A rule of 

thumb is that finite volume effects are small if L^m^ > 
4 — 5. For all quark masses used in this work we have 
LsTO^ > 4.6, except for dynamical Wilson fermions at the 
smallest quark mass where we have Ls^u- = 3.6. Since 
we do not have dynamical Wilson configurations on a 
larger volume we test for finite size effects in the hybrid 
scheme for which, at the smallest quark mass, there are 
MILC configurations for Ls — 2.5 and Lg = 3.5 giving 
LsTOtt = 4.6 and L^m^ = 6.4, respectively. In Fig. [9] we 
show results for the N to A axial form factors G^[Q^) 
and C|f (Q^) as well as GttAtaCQ^) for these two spatial 
sizes. Results on the smaller lattice are consistent with 
results on the larger lattice. This indeed shows that finite 
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FIG. 9: The upper graph shows C^, the middle graph 
and the lower graph GnNA as a function of in the hybrid 
approach for spatial volumes, 20^ (crosses) and 28"^ (open 
triangles) at the smallest quark mass, mi — 0.01. 



volume effects are small for Lsm.^ > 4.5. Since for all our 
quark masses, except the lightest mass dynamical Wilson 
fermions LsTtIt^ > 4.6, we expect finite volume effects to 
be small. We note, however, that a systematic study 
of volume effects that would allow an extrapolation of 
our quantities to infinite volume requires results using at 
least three different volumes. This is beyond the scope 
of the present work. 

Finally we comment on the evaluation of the kincmat- 
ical factors in the expressions for the form factors, which 
involve the masses of the nucleon and A and their en- 
ergies. The masses are evaluated using two-point func- 
tions in the standard way. The energies are calculated 
using the continuum dispersion relation, E = rrp- -I- p^, 
where — ni^'nlLs)'^^ n = 1,2,... One can com- 
pare results obtained using continuum dispersion rela- 
tions to those obtained with the lattice dispersion rela- 
tion sinh^(i<^) = sinh^m + ^j^-^ g sin^ 27rn,;/is. We 
find that the mean value of the form factors is almost 
unchanged. The Q^-values are also very close for 
1 GeV^. At larger momentum transfers the lattice dis- 
persion relation shifts the to smaller values. Using 
two-point functions to extract the energy also yields con- 
sistent results for the form factors albeit with larger er- 
rors. In what follows we will present results as a function 
of calculated using the continuum dispersion relation. 



V. RESULTS 

We first discuss results on quantities and ratios for 
which the renormalized quark mass is not required. This 
eliminates one source of systematic error, namely lattice 
artifacts on the value of m^. Furthermore, in general, ra- 
tios show weaker dependence on quark mass. For these 
reasons, they are more suited for comparison with phys- 
ical results. 
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FIG. 10: The ratio of form factors G.^na{Q^) /G^nn{Q^) 
as a function of for Wilson fermions for the quenched 
theory, denoted by JVf = 0, at k = 0.1554 (m^ = 0.56 GeV), 
K = 0.1558(m^ = 0.49 GeV) and k = 0.1562(m^ = 0.41 GeV) 
and for two dynamical Wilson quarks, denoted by Nf = 2, 
at K = 0.1575 (m„ = 0.69 GeV) [H, k = 0.1580 {m„ = 0.51 
GeV) [23] and K = 0. 15825 (m^ = 0.38 GeV) [H. The dashed 
line is the result of fitting the quenched results to a constant, 
yielding a value of 1.60(2). 
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FIG. 11: The ratio of 2C|*(Q2)/Ga(Q^) as a function of Q^. 
The notation is the same as that of Fig. 1101 Fitting the 
quenched results to a constant yields a value of 1.63(1). 

For the same lattice momentum vectors the values 
for the nucleon system differ from those in the iV — A sys- 
tem. In order to take ratios of form factors computed in 
these two different systems we interpolate the form fac- 
tors in the nucleon system to the value of the N — A 
system. In Fig. [10] we show the ratio of the form factors 
GTrNA{Q'^)/GTrNN{Q^) for qucnchcd and two degenerate 
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FIG. 12: The ratio of SCi (Q^) /Gp{Q^) as a function of 
The notation is the same as that of Fig. 1101 Fitting the 
quenched results to a constant yields a value of 1.73(3). The 
dotted line denotes the value of 1.60 obtained by fitting the 
ratio Gttn i^iQ'^) / Gt,n n{Q'^) ■ 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



"T 1 1 1 1 r 



X N =0 m =0.56 GeV 

P IT 

O N =0 m =0.49 GeV 

F ir 

* N =0 m =0.41 GeV 



N,=2 m^=0.69 GeV 
• N =2 m =0.51 GeV 

F It 

..p=^ m^=U.li*j ue, 
I !_! 1_ 



0.2 



0.4 



0.6 0.8 1.0 

Q2 (GeV^) 



1.2 



1.4 



FIG. 13: The ratio of Eq. 112} to Eq. HI}. The notation is 
the same as that of Fig. 1101 



flavors of dynamical Wilson quarks denoted by Np — 
and Np = 2 respectively. As can be seen, this ratio 
is independent and shows no statistically significant 
quark mass dependence. Fitting the quenched results 
to a constant we obtain the value of 1.60(2) shown by 
the dashed line. If we assume pole dominance for the 
form factors Gp{Q^) and Cq{Q'^) then the GTRs sim- 
plify to the relations given in Eq. (fH|) . Taking the ratio 
of the diagonal and non-diagonal relations we find that 
G^na{Q^)/G^nn{Q^) = 2Ci[Q^)/GA{Q^)- In Vig.m 
we show the ratio 2C^ {Q^) / G a{Q^) , which is indeed also 
independent, and fitting to the quenched data we find 
the value of 1.63(1) shown by the dashed line. Therefore, 
on the level of ratios, the GTRs are satisfied. We can use 
the relations given in Eq. (|13p for Gp and Gq to eliminate 
Gttnn and G-^na in Eqs. (fTT|) and to obtain 



GpiQ') 



Am%/ml 



1 + Qyml 



Ga(Q') 



G^iQ' 



/ml 



1 + Q'^/ml 



Gi(Q^) 



(47) 



These relations are again a manifestation of pion pole 
dominance. Taking ratios, we find that %G^{Q'^) /Gp{Q'^) 
should be equal to the ratio 2C^{Q^)/Ga{Q^) and con- 
sequently to G7rjvA(Q^)/G^ArAr((5^)- As can be seen in 
Fig. [T^ we indeed find that also this ratio is constant as 
a function of Q^. Fitting the quenched data to a constant 
we obtain the value of 1.73(3), shown by the dashed line 
in the figure. This is about 6% larger than what we find 
for the other two ratios. Therefore we conclude that ra- 
tios based on the relations given in Eq. (fH|) are better 
satisfied than ratios obtained using Eq. (|T7)) . Relaxing 
the assumption on pion pole dominance of Gp and Gq, 
we can consider directly the ratio of the non-diagonal 
to the diagonal GTR given in Eqs. p^ and Eqs. pTj) re- 
spectively. As can be seen in Fig. [131 the ratio is indeed 
consistent with unity. 
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FIG. 14: The 
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ratio of nucleon axial form factors 
Gp{Q'') / G a{Q^) for Wilson fermions for the quenched theory 
and for two dynamical Wilson quarks using the same notation 
as in Fig [To] We also show results from Ref. [l^l obtained in 
the hybrid approach using the same quark masses as the ones 
used in this work, namely mi — 0.03 (star), mi — 0.02 (filled 
triangle) and mi = 0.01 (inscribed squares). The dash line 
shows the expected behavior assuming pion pole dominance 
as given in Eq. (|47[) . where for m-^ and mjv we use the values 
computed on the lattice at k = 0.1562. The solid curve is a 
fit to a monopole form of the quenched data at k = 0.1562. 

Having examined the Q^-dependence on the level of 
ratios of GTRs in the nucleon and N-A systems we now 
discuss the dependence of the form factors for the two 
systems separately. In Fig. [U we show the ratio of nu- 
cleon axial form factors Gp{Q^) /Ga{Q'^) as a function of 
for quenched and two degenerate fiavors of dynamical 
Wilson quarks. Recent results from Ref. [l^] obtained 
in the hybrid scheme at the same quark masses as the 
ones used in this work for the calculation of the TV to A 
form factors are also included. In all cases, the ratio de- 
creases with confirming the stronger Q'^-dependence 
expected for Gp{Q^) as compared to Ga{Q^)- A simi- 
lar behavior is also observed for the corresponding ratio 
C^(Q2)/C^(q2) N to A shown in Fig. M If pion 
pole dominance holds, then the ratios Gp{Q^)/Ga{Q^) 
and Gq{Q'^)/G^{Q^) should be described by the rela- 
tions given in Eq. ([47]) with no adjustable parameters. 
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FIG. 15: The ratio of to A axial transition form factors 
Ce{Q^)/CsiQ'^)- The notation is the same as that of Fig. ITU 
The dotted line shows the prediction of pion pole dominance 
predicted in Eq. (|47p but for the hybrid case at the lightest 
quark mass. 
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FIG. 16: The upper graph shows Ga{Q^) and the lower graph 
Gp{Q'^) as a function of Q^. The solid curve is a fit to a 
dipole form of the quenched results at «: = 0.1562. The fit to 
an exponential form shown by the dashed line falls on top. 
The dotted line shown in the upper graph corresponds to a 
dipole form with axial mass tua — 1.1 GeV used to described 
experimental data. The dashed line in th elower graph shows 
the result expected from pion pole dominance in Eq. (|47|l . 
The solid line corresponds to Eq. 1501 The rest of the notation 
is the same as that in Fig. 1141 Results in the hybrid approach 
are from Ref. [13]. 

I 



In order to examine the Q^-dependence of the form 
factors separately and compare with contimium quan- 



In Figs. [13] and [15] we show, with the dashed lines, the 
resulting curves for the case of quenched lattice results 
at the lightest pion mass obtained assuming the relations 
given in Eqs. l[47|) . As can be seen, as 0, both ra- 

tios increase slower than pion pole dominance predicts. 

In addition we show by the solid curves fits to the same 
1.0| ^ , , 1 
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FIG. 17: The upper graph shows C^{Q^) and lower graph 
C^(Q^) as a function of . The notation is the same as that 
in Fig. [H 



quenched data using a dipole form 



(g2/m2 + l) ' ' 

with fit parameters cq and m. The values of co and ni 
extracted from the fits are given in Table III. In the 
quenched case, we find that m > m^r, whereas for the 
hybrid scheme, although m ^ 771^^, cq is smaller than 
m%/m1 causing the dotted Hne shown in Fig. [TH] ob- 
tained at the lightest quark mass, to be higher than the 
corresponding lattice results. The calculation in the hy- 
brid scheme at the lightest quark mass is done on a larger 
lattice enabling us to compute the form factors at low 
Q^-values, much lower than in the case of dynamical Wil- 
son fermions. These results show clear deviations from 
quenched results at low Q^, where pion clouds effects are 
expected to dominate. 



I 

titles, we need to multiply lattice results with the axial 
renormalization constant Za ■ These constants are known 
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Nucleon elastic 


(GeV) 


m (GeV) 


Co 


niA (GeV) go tHa (GeV) go A 


givNN 


Quenched Wilson fermions 


0.563(4) 
0.490(4) 
0.411(4) 


0.671(14) 
0.597(14) 
0.511(16) 


13.71(34) 
15.23(43) 
17.70(76) 


1.659(20) 1.088(8) 1.271(9) 1.074(5) 0.110(2) 
1.632(19) 1.079(7) 1.249(9) 1.069(5) 0.083(2) 
1.578(28) 1.080(12) 1.220(10) 1.066(6) 0.062(2) 


9.943(99) 
9.126(93) 
8.410(100) 


10.609(73) 
10.143(91) 
9.725(140) 


Nf ~ 2 dynamical Wilson fermions 


0.691(8) 
0.509(8) 
0.384(8) 


0.750(43) 
0.642(77) 


14.13(1.01) 
11.15(1.82) 


1.831(22) 1.067(6) 1.393(16) 1.063(6) 0.114(3) 
1.709(46) 0.999(17) 1.296(29) 0.995(17) 0.038(15) 
2.019(78) 0.951(18) 1.528(44) 0.943(15) 0.044(10) 


11.48(245) 


10.486(122) 
9.071(294) 
8.613(551) 


Nucleon to A 


(GeV) 


m (GeV) 


Co 


ruA (GeV) 50 ifiA (GeV) go A' 






Wilson fermions quenched 


0.563(4) 
0.490(4) 
0.411(4) 


0.691(13) 
0.624(15) 
0.545(16) 


3.44(79) 
3.75(12) 
4.23(17) 


1.544(32) 0.952(16) 1.205(8) 0.926(5) 0.106(2) 
1.537(33) 0.930(14) 1.192(10) 0.910(7) 0.079(2) 
1.534(36) 0.906(15) 1.189(13) 0.887(9) 0.052(2) 


16.560(194) 
14.692(188) 
12.609(180) 


17.174(166) 
16.195(206) 
14.873(264) 


Wilson fermions, dynamical Np = 2 


0.691(8) 
0.509(8) 
0.384(8) 


0.604(95) 
0.352(151) 
0.379(58) 


4.75(1.04) 
8.38(6.11) 
6.34(1.69) 


1.696(51) 0.988(24) 1.368(13) 0.937(5) 0.109(2) 
1.760(59) 0.865(25) 1.454(46) 0.808(20) 0.063(2) 
1.968(118) 0.843(40) 1.410(51) 0.808(24) 0.024(15) 




17.536(190) 
14.970(452) 
12.685(1.416) 


0.594(1) 
0.498(3) 
0.357(2) 


0.576(28) 
0.485(27) 
0.389(18) 


5.08(22) 
6.15(52) 
8.59(64) 


Hybrid action 
1.924(85) 0.883(22) 1.477(42) 0.868(15) 0.076(5) 
1.892(101) 0.864(32) 1.505(71) 0.835(27) 0.0648(7) 
1.849(71) 0.760(18) 1.522(72) 0.708(19) 0.054(5) 


10.053(324) 


17.649(236) 
17.329(496) 
14.815(558) 



TABLE III: The first column gives the pion mass in GeV, the second and third columns the fit parameters m and cq extracted 
from fitting the ratio Gp/Ga {C^/C^) for the nucleon (N to A) case, the fourth and fifth columns the dipole parameters 
tua and go extracted from fitting Ga {C^) for the nucleon (N to A) and the sixth and seventh columns the corresponding 
parameters but using an exponential Ansatz go exp( — Q'^/m^). The eight column gives A or A' defined in Eq. The last 

two columns give the value of the strong coupling constants gTvNN or g-nNA- The first value of the strong coupling constant is 
determined using the fit function of Eq. (f55)l , whereas the second using a linear fit according to Eq. (f5^ . 



for both Wilson fermions and DWF within the hybrid 
scheme. The values that we use are given in Table IV. 
We collect our lattice results for the nucleon form factors 
in Tables V and VI and for to A in Tables VII, VIII 
and IX of the Appendix. All errors are calculated using 
jackknife analysis. In Fig. [111 we show Ga{Q^) and the 
induced pseudoscalar form factor Gp{Q^). For compari- 
son we also show results obtained in the hybrid approach 
from Ref. [l^]. The main observation is again that at 
the smallest domain wall quark mass, the hybrid results 
show deviations. In particular, we note that the value 
for the nucleon axial charge gA becomes larger in the hy- 
brid scheme approaching the experimental value. This is 
in agreement with the findings of Ref. j^] . Since there 
are recent state-of-the-art lattice studies of gA [H, HSl 
we will not discuss it further here but rather investigate 
the dependence of the form factors. We also find that 
Gp(Q^) increases more rapidly at low in the hybrid 
scheme when the pion mass decreases to about 350 MeV. 
In Fig. [T7] we show the corresponding N to A transition 
form factors C^{Q^) and Cq{Q'^). The hybrid results 
show the same behavior as in the case of the nucleon form 



factors, yielding a different behavior at low when the 
pion mass becomes about 350 MeV. The Q^-dependence 
of both GaIQ^) and G^{Q^) can be well described by a 
dipole Ansatz 



This is what is usually used to described experimental 
data for Ga{Q^) where a value of niA ^1.1 GeV is ex- 
tracted for the axial mass. The same dipole Ansatz is 
also used to described C^(Q^), where an axial mass of 
1.28±0.10 GeV [Si has been found. In addition, we fit to 
an exponential form given by goe~'^ . Both Ansatze 
describe well our results as can be seen in Figs. [16] and 
[T7I where the two lines, which are fits to quenched lat- 
tice results at the smallest quark mass, can hardly be 
distinguished. The values of the axial masses extracted 
from these fits are given in Table III. We find an axial 
mass that is larger than what is deduced from exper- 
iment. This means that Ga{Q'^) and C^{Q'^) fall off 
slower than in experiment. This is clearly seen in Fig. 1161 
where we include the dipole curve taking uia = 1.1 GeV. 
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K or mi 



afn/ZA 



Za 



Quenched Wilson fermions 



0.1554 0.0403(4) 0.0611(14) 0.808(7) [42| 

0.1558 0.0307(4) 0.0587(16) 0.808(7) 

0.1562 0.0213(4) 0.0563(17) 0.808(7) 

Nf = 2 Wilson fermions 

0.1575 0.0441(4) 0.0649(8) 0.77(2) [43] 

0.1580 0.0229(4) 0.0494(9) 0.78(4) [43] 

0.15825 0.0122(3) 0.0467(13) 0.8° 



Hybrid action 
0.03 0.0475(3) 0.0678(6) 
0.02 0.0324(4) 0.0648(8) 
0.01 0.0159(2) 0.0636(6) 



1.1085(5) [23] 
1.0994(4) [23] 
1.0847(6) [23] 



"Estimated from the values of Za at k = 0.1575 and 0.1580 

TABLE IV: The first column gives the hopping parameter k 
for Wilson fermions or the mass of the domain wall fermion, 
the second the renormalized quark mass, the third the un- 
renormalized pion decay constant /t^/Za in lattice units, and 
the fourth the axial renormalization constant Za- 



Having fitted Ga{Q^) and C^{Q'^), the Q^-dependence 
for the form factors Gp{Q^) and Cq{Q^) can be obtained 
using Eq. (|47p . The resulting curves are shown by the 
dashed line in Figs. [H] and [T7] and show deviations at 
low Q^. In addition we show curves that correspond to 



.90 Co 



(g2/mji + l)2(Q2/m2 + l) 



(50) 



with m extracted from fitting the ratio of 
Gp{Q^)/Ga{Q^) in the case of the nucleon system 
and C^(Q2)/c^(q2) foj. ^i^g expected this 

provides a good description of the (52_(jependence for 
both Gp{Q'^) and C^iQ"^) shown by the soHd Hues, 
which correspond to the parameters of the quenched 
data at K = 0.1562. 

We now present results that require knowledge of the 
renormalized quark mass. The renormalized quark mass, 
rUq, is determined by evaluating the pion to vacuum ma- 
trix element of the axial Ward-Takahashi identity given in 
Eq. ([5]). As mentioned already, for Wilson fermions the 
axial Ward identity is satisfied only up to 0{a) terms. 
We expect these corrections to become more severe as 
we approach the chiral limit. As we already mentioned, 
the quark mass in the hybrid scheme, , was tuned 
to reproduce the mass of the lightest pion in the stag- 
gered theory. Given that domain wall fermions satisfy 
the AWI when the size of the fifth dimension is taken 
to infinity, corrections to Eq. ([8|) come from the residual 
mass due the finite length of the fifth direction. There- 
fore differences between the values of and are 
due to chiral symmetry breaking becuase of the finite size 
of the fifth dimension. We show in Fig. [181 the renormal- 
ized quark mass extracted from the axial Ward identity 
for quenched and two dynamical Wilson fermions and in 
the hybrid approach. As can be seen, the pion mass ex- 
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FIG. 18: The renormalized quark mass iriq versus m„ for the 
quenched theory (crosses) , for two dynamical Wilson fermions 
(open circles) and for the hybrid scheme (filled squares). In 
the hybrid case we also show (asterisks) determined by 

tuning the pion mass [33]. The lines are linear fits to m^. 
The open square shows the extrapolated value of in the 
hybrid scheme. 
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FIG. 19: The upper graph shows the ratio 
/7rG,r]VJv((5^)/mjvG'A(Q^) and the lower graph the ra- 
tio /7rG^jvA(Q^)/2mjvC^(Q^). The notation is the same as 
in Fig. [ft] 



trapolates to zero at = for the quenched theory. 
For dynamical Wilson fermions m^r is not exactly zero 
at TUq — indicating finite a-corrections. For the hybrid 
scheme we show both the renormalized mass computed 
using the AWI, m^, and . As can be seen m^, is ap- 



proximately equal to confirming that the residual 

mass is small. The biggest deviation, as expected, is ob- 
serve at the smallest value of the quark mass. From these 
results we confirm that in the hybrid scheme oc niq 
within our statistical errors. 

Having determined ruq and /tt using Eq. ([40]) 
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FIG. 20: The upper graph shows the ratio 
AmNf7rG*ffff(Q^)/m^Gp{Q'^) and the lower graph the 
ratio niN fnG*N^{Q^)/2m,^Ce{Q^) for the hghtest quark 
mass considered in each of our three types of calculations. 
We have defined GInn{Q^) = G^nn{Q'^)/{1 + Q'^ /ml) with 
a corresponding expression for G^jvA- 



we can evaluate the form factors G-nNNiQ"^) and 
GttNa{Q^)- We first examine the Goldberger- 
Treiman relations as given in Eq. (|14p by con- 
sidering the ratios f-KGT^NN{Q^)lmMGA{Q^) and 
/^G^^^((5^)/2mArC^((5^), which should be equal to 
unity. Note that in these ratios the axial and pseu- 
doscalar renormalization constants cancel. These ratios 
are shown in Fig. [121 Wc find that, in the quenched 
theory, they are less than one for small but become 
one for ~ 0.5 GeV^. This is also approximately true 
for dynamical Wilson for the two heaviest quark mass. 
Results in the hybrid approach, on the other hand, show 
smaller deviations from unity at low . We also expect 
that the ratios 



m2Gp(g2) (l + QVm2) 



2mlCi{Q^){l + Qyml 



(51) 

should be unity if pion pole dominance is valid. As can 
be seen in Fig. 1201 for the lightest quark masses in our 
three types of action these ratios are indeed consistent 
with one. 

Finally we discuss the Q^-dependence of the form fac- 
tors GyrNwiQ'^) and CttAtaIQ^) separately. In Fig. [HH 
we show GTrNNiQ^) and Gt^na{Q^) at the smallest 
quark mass in the quenched theory and for dynami- 
cal Wilson quarks. For G.,rM\{Q^) we also show re- 
sults in the hybrid scheme at a similar quark mass. 
Results for these form factors at small are consis- 
tent with each other unlike the other form factors, indi- 
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FIG. 21: The upper graph shows Gttnn{Q^) for Wilson 
fermions at the smallest pion mass. The lower graph shows 
Gttna{Q^) for Wilson fermions and DWF at the smallest 
pion mass. The dashed lines follow from the GTRs relations 
given in Eq. (|14|l . The solid lines are fits using Eq. ((52). In 
the case of G',r]VA(Q^), we also show by the dashed-dotted 
line (with larger slope) the curve that corresponds to taking 
G-nN^iQ^) = 1.6 G„nn{Q^)- 



eating that unquenching effects on these quantities are 
small. Assuming PCAC and pion pole dominance, the 
Q^-dependence of Gt;nn{Q'^) and G-j^NAiQ^) is com- 
pletely determined from the GTRs given in Eq. (fH)) 
once we know Ga{Q^) and C^{Q^). Using the dipole 
Ansatz of Eq. ([i^ for the Q^-dependence of Ga{Q'^) and 
C^(Q^) with the parameters given in Table III, we ob- 
tain the dashed lines shown in Fig. [^Ij The discrepancy 
already observed in the ratio f7rG.,rNNiQ^)/iTi]yGA{Q^) 
and /,rG^ArA(Q^)/wArC^((5^) at low values is clearly 
seen here. Results in the hybrid approach confirm devi- 
ations from the GTRs at low Q^. The Q^-dependence 
of GyrNNiQ"^) and GjrNAiQ^) can be described using a 
linear Ansatz given by 



GttNa{Q^ 



a'(l-A'^ 
mi 



(52) 



with a (a') and A (A') fit parameters. These linear fits 
are shown by the solid curves in Fig. [^T] and provide 
a good description to the results. Note that we have 
excluded from the fits the value at the lowest in all 
cases in the quenched theory as well as in the hybrid 
approach at the smallest quark mass, since our statistical 
error maybe underestimated for this value of Q^. The 
values we find for the parameters a (a'), which determine 
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FIG. 22: The upper graph shows G^NNiQ'^) and the lower 
graph shows GnNA{Q^) as a function of . The dashed hnes 
are fits to the quenched results at k = 0.1562 obtained using 
the functions of Eq.[50]and allowing an overall constant to be 
fitted. 



the strong coupling constant g^^NN {9tiNa), and A (A') 
are given in Table III. Note that A and A' decrease as 
the quark mass decreases. In the quenched theory at 
the smallest quark mass A ~ 6%. One expects that 
this value decreases further as we approach the physical 
limit becoming comparable to the value of A = 2.44% 
obtained using baryon chiral perturbation theory p^ . 
However, the corresponding value of a is smaller than 
the expected value of mj^gA/ f-n and decreases with the 
quark mass being about 83% less at the heaviest and 
about 73% at the lightest quark mass in the quenched 
theory. The relation Gt^na[Q^) — l-6G'7rAfjv(Q^) can be 
used to determine the behavior of Gt^na from that 
of Gttnn ■ The dashed-dotted line in Fig. [21] shows the 
resulting curve, which approximates the results well. 

In Fig. mi we show results for G-j^NAiQ'^) and 
Gttnn{Q^) at all values of the quark masses considered 
in this work. Assuming pion pole dominance we relate 
the Gttnn and G-j^na to Gp and Gg respectively using 
Eq. (gT]). Taking the functions defined in Eq. for the 
Q^-dependence of Gp{Q^) and Gq{Q'^) , we write 



G-r^NNiQ^) = K 



N 



Q^/ml + 1 



(Q2/TO^ + l)2(Q2/m2 + l) 



(53) 



with a corresponding expression for G-^na (Q^ ) • The only 
free parameter is an overall constant, Kn, to be fitted to 
the results. The value of Kn determines Qtinn and is 
given in Table III. The fits to the quenched data at the 
smallest quark mass using Eq. ([55)) are shown by the 
dashed lines in Fig. [22| Note that if m were the pion 
mass, then the pole term would cancel, leaving a dipole 



(52-dependence for these form factors. Allowing m ^ TO,r 
and adjusting the overall strength, we can obtain a rea- 
sonable description of the Q^- dependence of CttAt n and 
Gttna- Since lattice results for these form factors do 
not increase at low values of Q'^ as fast as expeced by 
PCAC, we obtain a smaller value at = 0. The values 
of 9TrNN ~ G^Ar7v(0) and g-rrNA — G-^na{0) extracted 
using Eq. (|53p. when the fit yields x^/d.o.i. ^ 1.5, are 
given in Table III. If one uses the relation G-j^NNiQ'^) = 
/ f-R)G AiQ"^) as the GTR would suggest, then the 
extrapolated value at the lightest pion mass would be 
girNN — 11-8 ± 0.3 in the quenched theory closer to the 
experimental value of 13.2i;to'o5 L45J- Therefore the dif- 
ferent low Q2 dependence observed in the lattice results 
compared to what is usually assumed, is responsible for 
the lower values of Qt^nn and Qt^na extracted from these 
fits. 



VI. SUMMARY AND CONCLUSIONS 

We have presented results for the nucleon axial-vector 
form factors Ga{Q^) and Gp{Q^) as well as for the 
corresponding N to A axial transition form factors 
C^(Q2) and CiiQ^). The nNN and nNA form fac- 
tors Gt;nn{Q'^) and G-^NAiQ^) are also evaluated. Us- 
ing ratios that show very weak quark mass dependence 
and are therefore expected to have the same value in the 
physical limit, we reach a number of phenomenologically 
important conclusions. One of the main conclusions is 
that GttNN and Gt^na have the same dependence 
yielding a ratio of G-^NAiQ"^) /GTrNNiQ^) — 1-60(2) in 
good agreement with what is expected phenomenologi- 
cally. Similarly the ratio 2C^(Q2)/g^(q2) ^ 1.63(1) 
is also independent of Q2 Equality of these two ra- 
tios implies the Goldberger-Treiman relations. The ra- 
tio 8Gq{Q'^)/Gp{Q'^) on the other hand is larger by 
about 6%. The popular pion pole dominance hypoth- 
esis is examined using our lattice results. We find that in 
the quenched theory the ratios of Gp{Q'^)/Ga{Q'^) and 
Gg{Q'^)/G^{Q'^) require a larger pole mass m than the 
corresponding pion mass in order to get a good descrip- 
tion at low values of . On the other hand, these ratios 
in the hybrid approach are well described with m m^. 
However the overall strength differs from what is pre- 
dicted using the Goldberger-Treiman relations. 

We also studied the Q2_(igpendence of the form fac- 
tors separately. Comparing quenched and unquenched 
results at pion mass of about 350 MeV, we observe large 
unquenching effects on the low Q2_(igpg]^(jence of the four 
form factors, Ga(Q^), Gp{Q^), C^{Q^) and Gi{Q^). 
This confirms the expectation that pion cloud effects 
are expected to be large at low Q^. We find that the 
(52-dependence of the form factors Ga and is well 
described by a dipole of the form co/(l + Q'^/m\)'^. 
For the pion masses used in this work we find that 

niA ^ 1-5 GeV as compared to 1.1 GeV used to de- 
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scribe the of experimental data for Ga- An expo- 
nential Ansatz of the form cq exp(— Q^/m^) also pro- 
vides a good description of the Q^-dependence. In 
agreement with a recent lattice evaluation of the nu- 
cleon axial charge [23 |, we find that qa increases 
when using unquenched configurations on a large vol- 
ume and pion mass of about 350 MeV becoming con- 
sistent with experiment. A different low Q^-dependence 
is observed for both quenched and unquenched results 
in the case of G-j^NNiQ"^) and Gt^na{Q^)- Instead of 
the dipole form expected from the Goldberger-Treiman 
relations of Eq. (fT4|) one finds that Gt^nn and G^ata, 
in the limit ^ increase less rapidly. As a result 
of this, the values that we extract for the strong cou- 
pling constants g^^NN = limQ2^o G^NNiQ'^) and g-^NA = 
limg2^Q GTrNAiO) are smaller than those extracted from 
experiment. One ingredient that is needed for the deter- 
mination of these form factors is the renormalized quark 
mass. In this work, we use the axial Ward identity to 
determine it. On the lattice, the axial Ward identity has 
0{a) corrections in the case of Wilson fermions, which 
can become important in particular as we approach the 
physical limit. In the case of domain wall fermions, the 
axial Ward identity is only exact in the limit of large 
fifth dimension with corrections due to the residual mass, 
which again become more important in the chiral limit. 



The renormalized quark mass, however, affects the over- 
all strength of these form factors and therefore it can not 
explain the different Q^-dependence. To investigate this 
issue further one would like to use lighter quark masses 
on a finer lattice. This will become feasible in the near 
future as such dynamical simulations are under way. 
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VII. APPENDIX 



Nucleoli elastic:Quenched Wilson fermions 



(OcV^) 


GaIZa 


Gp/Za 


GttNN 




in^ = o.r,b 


(Gc\-) 




0.0 


1.332(17) 






0.17 


1.193(13) 


12.117(188) 9.645(232) 


0.34 


1.077(16) 


8.456(169) 


9.296(139) 


0.49 


0.977(15) 


6.227(153) 


8.804(202) 


0.64 


0.893(15) 


5.102(157) 


8.294(240) 


0.79 


0.825(18) 


4.123(95) 


7.799(254) 


0.93 


0.752(18) 


3.277(111) 


7.318(286) 


1.19 


0.662(25) 


2.502(117) 


6.352(374) 


1.32 


0.589(27) 


2.037(102) 


5.738(380) 


1.44 


0.557(31) 


1.885(114) 


5.384(456) 


1.56 


0.502(32) 


1.563(106) 


4.910(475) 


1.68 


0.421(38) 


1.220(141) 


3.780(565) 


1.79 


0.451(47) 


1.307(147) 


4.177(590) 


1.90 


0.371(40) 


1.001(119) 


3.374(497) 


2.12 


0.337(86) 


0.951(274) 


2.437(881) 


m,r = 0.490(4) (GeV) 


0.0 


1.325(15) 






0.17 


1.183(15) 


12.432(321) 


9.109(190) 


0.33 


1.065(21) 


8.352(147) 


8.845(166) 


0.49 


0.961(23) 


5.964(169) 


8.419(243) 


0.64 


0.879(20) 


4.914(166) 


7.850(247) 


0.78 


0.811(15) 


3.902(116) 


7.457(262) 


0.91 


0.738(22) 


3.044(109) 


7.065(318) 


1.17 


0.654(28) 


2.317(124) 


6.183(421) 


1.29 


0.577(30) 


1.874(108) 


5.534(428) 


1.41 


0.546(33) 


1.750(119) 


5.037(487) 


1.53 


0.489(36) 


1.426(110) 


4.677(518) 


1.64 


0.409(42) 


1.094(147) 


3.644(657) 


1.75 


0.451(55) 


1.234(161) 


4.055(651) 


1.86 


0.363(45) 


0.917(126) 


3.246(537) 


2.07 


0.323(93) 


0.853(282) 


2.033(903) 


m,r = 0.411(4) (GeV) 


0.0 


1.319(19) 






0.17 


1.177(19) 


12.824(345) 


8.619(182) 


0.33 


1.054(22) 


8.225(202) 


8.417(232) 


0.48 


0.947(19) 


5.649(185) 


8.085(277) 


0.63 


0.867(28) 


4.700(199) 


7.380(309) 


0.76 


0.798(21) 


3.650(117) 


7.102(301) 


0.90 


0.723(27) 


2.762(120) 


6.886(393) 


1.14 


0.647(36) 


2.084(142) 


6.167(535) 


1.26 


0.566(35) 


1.685(117) 


5.405(550) 


1.38 


0.538(40) 


1.601(128) 


4.601(541) 


1.49 


0.474(41) 


1.269(118) 


4.434(589) 


1.60 


0.394(47) 


0.938(164) 


3.549(854) 


1.70 


0.453(67) 


1.148(187) 


4.009(777) 


1.81 


0.353(52) 


0.817(142) 


3.217(633) 


2.00 


0.324(115) 


0.778(322) 


1.578(986) 



TABLE V: The first column gives the in GcV^, the second 
Ga, the third Gp and the fourth G^nn for quenched Wilson 
fermions. The errors quoted are jackknife errors. 
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Nucleon elastic: Nf — 2 Wilson fermions 



(GeV^) 


Ga/Za 


Gp/Za 


G-^NN 


= 0.691(8) (GeV) 


0.0 


1.382(11) 






0.43 


1.105(16) 


8.976(298) 


9.367(228) 


0.82 


0.893(16) 


5.115(170) 


8.361(278) 


1.19 


0.735(22) 


3.075(172) 


7.696(361) 


1.53 


0.658(36) 


2.575(193) 


6.621(516) 


1.86 


0.592(44) 


2.033(190) 


5.980(552) 


2.17 


0.484(45) 


1.325(159) 


5.460(617) 


2.75 


0.356(91) 


1.108(322) 


2.441(788) 


3.03 


0.290(74) 


0.678(212) 


3.081(853) 


= 0.509(8) (GeV) 


0.0 


1.270(27) 






0.42 


0.982(21) 


7.229(388) 


8.505(350) 


0.81 


0.809(24) 


3.914(221) 


7.827(331) 


1.16 


0.657(32) 


1.983(225) 


7.913(517) 


1.49 


0.502(50) 


0.668(218) 


7.689(1.057) 


1.79 


0.457(43) 


0.663(151) 


7.005(860) 


2.09 


0.377(72) 


0.300(161) 


7.416(1.691) 


2.63 


0.216(74) 


0.702(143) 


4.551(1.841) 


2.88 


0.104(60) 


-0.025(84) 


3.298(1.897) 


3.12 


0.028(25) 


0.226(54) 


0.651(0.442) 


= 0.384(8) (GeV) 


0.0 


1.175(33) 






0.42 


0.990(25) 


5.418(370) 


7.812(391) 


0.78 


0.835(32) 


3.395(244) 


6.454(366) 


1.11 


0.786(49) 


2.269(221) 


5.166(549) 


1.41 


0.611(50) 


1.270(220) 


5.673(689) 


1.68 


0.632(98) 


1.361(223) 


4.867(1.017) 


1.95 


0.476(77) 


0.979(177) 


2.058(734) 


2.43 


0.092(527) 0.241(1.377) 0.638(5.365) 



TABLE VI; The first column gives the in GeV^ the second 
Ga, the third Gp and the fourth GttNN for Wilson fermions. 
The errors quoted are jackknife errors. 
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N to A: Quenched Wilson fermions 
(GeV^) C^IZa Cj/ZA G^na 
rrin = 0.563(4) (GeV) 



0.16 


1.016(14) 


2.675(81) 


15.613(356) 


0.35 


0.902(8) 


1.778(52) 


15.207(298) 


0.53 


0.803 


15) 


1.269(45) 


14.441(361) 


0.70 


0.712 


14) 


1.012(38) 


12.795(501) 


0.87 


0.635 


17) 


0.774(26) 


11.882(415) 


1.03 


0.571 


18) 


0.607(26) 


11.203(442) 


1.34 


0.465 


21) 


0.418(24) 


9.610(538) 


1.48 


0.414 


23) 


0.348(23) 


8.440(515) 


1.63 


0.371 


25) 


0.305(26) 


7.318(681) 


1.77 


0.330 


25) 


0.247(23) 


6.895(639) 


1.90 


0.295 


31) 


0.204(26) 


5.744(764) 


2.03 


0.273 


32) 


0.188(25) 


6.216(856) 


2.16 


0.230 


28) 


0.145(20) 


5.106(669) 


2.42 


0.172 


48) 


0.116(39) 


3.603(1.527) 




= 


0.490(4) (GeV) 


0.16 


0.999 


14) 


2.741(105) 


14.396(351) 


0.35 


0.885 


16) 


1.741(52) 


14.343(325) 


0.53 


0.787 


16) 


1.202(46) 


13.758(354) 


0.70 


0.691 


18) 


0.950(44) 


11.918(535) 


0.87 


0.616 


19) 


0.712(30) 


11.200(441) 


1.03 


0.553 


22) 


0.552(26) 


10.717(456) 


1.34 


0.450 


26) 


0.377(25) 


9.412(592) 


1.48 


0.401 


26) 


0.313(24) 


8.129(572) 


1.62 


0.357 


28) 


0.274(27) 


6.807(751) 


1.76 


0.315 


29) 


0.220(24) 


6.563(708) 


1.90 


0.288 


36) 


0.183(27) 


5.584(883) 


2.03 


0.262 


36) 


0.169(26) 


6.176(966) 


2.15 


0.219 


31) 


0.129(21) 


5.015(735) 


2.40 


U.ioU 


48) 


0.092(37) 


3.466(1.749) 




= 


0.411(4) (GeV) 


0.15 


0.975 


19) 


2.804(129) 


12.928(392) 


0.34 


0.864 


16) 


1.688(71) 


13.263(315) 


0.53 


0.769 


19) 


1.114(58) 


13.023(474) 


0.71 


0.668 


19) 


0.876(52) 


10.784(719) 


0.87 


0.593 


23) 


0.634(34) 


10.453(479) 


1.04 


0.536 


25) 


0.488(27) 


10.225(563) 


1.34 


0.438 


30) 


0.333(28) 


9.496(718) 


1.49 


0.391 


31) 


0.272(25) 


7.978(684) 


1.63 


0.347 


32) 


0.240(30) 


6.303(901) 


1.77 


0.303 


34) 


0.190(26) 


6.269(848) 


1.90 


0.284 


44) 


0.157(29) 


5.555(1.089) 


2.03 


0.250 


44) 


0.150(28) 


6.346(1.179) 


2.15 


0.211 


37) 


0.113(22) 


5.127(870) 


2.40 


0.121 


50) 


0.065(37) 


3.425(2.354) 



TABLE VII: The first column gives the in GeV^ the sec- 
ond C^, the third Cq and the fourth Gttna- The errors 
quoted are jackknife errors. 



N to A: Nf = 2 Wilson fermions 



(GeV^) 


C^/Za 


Cq/Za 


GttNA 


m,r = 0.691(8) (GeV) 


n AAQ 




z.uzzi^yo J 


lD.UoZl^40 / } 


U.ooZ 


u. / Do^^iy 


l.\Jo\JyO ( J 


lO.OOOl^4oU ) 


l.ZoO 


n finE;/'9Q\ 




IZ.ZDo^^DDO ^ 


i.DOD 




U.40 ( yAiO } 


1 n (^71 (QA ^^ 

lU.D ( 1 l^y40 ) 






n "^71 f 'iQ^ 
u.o ( i i^oy J 


o.yyui^ooo ) 


9 ^1^7 
Z.OO 1 




n 0Ar\f 
u.z^ui^oz J 


S 91 090^ 

o.ziui^i.uzy J 


o.UUo 


U.zzy^Oo^ 


U.l0z(4y J 


4.4oi(i.o41 J 


o.oiU 




n nsQ/' Q7^ 


4.yon^i.oyi ) 






tjuyio) ^^vjrt; V 








1 c^l 7^ 1 1 ^ 

1.01 t yiLL) 


lO. / OH^DZZ ^ 


u.oyo 


U.U 1 oi^zo 


n 7AA^^^^ 


ll.Dl 1 l^DDo J 


l.Oii 


u.oozi^zy^ 


U.OoUl^4D J 


y.yooi^i.U4y ) 


1 f{Q7 


u.ouzi^oo 


1 U.Z4Ul^0OJ 


o.oyyi^i.ouu j 


9 HPVQ 

z.uoy 


n ^ft9/'^8'i 

U.OOZ^^OO^ 


U.lOl l^Zo j 


7 ^SQ/'l 1 94^ 

1 .ooyi^i.iz4 J 


9 /ini 






7 91 Q/'l QQ^^ 


3.035 


0.144(62) 


0.039(25) 


3.656(2.079) 




ii'tt = 0. 


'J81(8) (Gc\') 




0.443 


0.837(28) 


1.335(96) 


11.302(722) 


0.891 


0.700(25) 


0.583(46) 


11.997(792) 


1.293 


0.620(34) 


0.381(35) 


8.956(1.0358) 


1.661 


0.530(38) 


0.256(29) 


9.168(1.482) 


2.002 


0.452(63) 


0.203(30) 


5.966(1.504) 


2.322 


0.342(48 


0.136(22) 


2.054(1.086) 


:!.i8;! 


0.i()2(2;'J) 


0.();U)(!)) 


().!):)()(().!)()()) 



TABLE VIII: The first column gives the in GeV^, 
second C^, the third and the fourth G^na- The ei 
quoted are jackknife errors. 
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N to A: Hybrid action 



(GeV^) 


C^/Za 


Cq /Za 


GttATA 


= 0.594(1) (GeV) 




(J.Doy(4^ J 


0101/^ AO\ 

z.iji^y^ ) 


lD.DUz(449j 




A /I 0/ OAA 

U.D4z(oUj 


l.ozDt^bUj 


lD.O09^06o) 


0.738 


0.578(28) 


0.944(57) 


14.796(808) 


0.983 


A /t AA /' 1 \ 

U. 490(31) 


0.614(43) 


13.944(1.049) 


1 91 


n 444/' '^n'l 


0.477(32) 


12.911(1.094) 




n ^Q9/''^9'l 

u.oyzi^oz J 


0.377(31) 


11.178(1.315) 


1.0 1 41 


U.Ol / l^'iZ J 


0.244(45) 


10.628(2.508) 




n 9fii /'/1 1 

U.Zult 1 


0.180(32) 


9.296(2.253) 


9 97Si 
Z.Z ( 


n 9nn/'47'i 


0.132(29) 


9.768(3.075) 


9 A79 


n 901 f 


0.118(32) 


10.702(3.698) 


9 RRH 
Z.DOU 


U. iOZl^oU } 


0.115(63) 


4.125(4.734) 


9 

Z.O'y:'^: 


Dim f 9i.Q\ 


0.041(47) 


11.147(10.192) 




A 


498(3) (GeV) 


A 1 ni 
U.iyi 


U.Doo^oo J 


2.398(186) 


15.404(608) 


0.471 


0.621(27) 


1.212(77) 


15.823(667) 


0.735 


A Cr/tA/'0'~l\ 

U.540(3ij 


0.809(75) 


15.325(896) 


u.yoo 


n c;4s/'49^ 


0.686(76) 


13.718(1.466) 


1 99/1 




0.445(44) 


12.040(1.226) 


1 zl^9 




0.332(40) 


10.760(1.523) 


1 SS9 


u.zyoi^oo J 


0.191(42) 


10.334(3.357) 


9 ns7 

Z.Uo ( 


A OQ/l /(^T^^ 

U.Z(54|^D ( J 


0.182(47) 


4.914(2.902) 


9 9^^4 


'^9'^f'l 1 i\ 
U.OZOl^ 1 lU 


0.190(70) 


2.688(5.406) 


9 /17R 


U.ZUol^lUD 


0.101(59) 


2.967(5.431) 


9 i^i^n 
z.oou 


n 1 9 Pi n ^4 


0.102(99) 






™, A 


357(2) (GeV) 


A ACA 


A Ci?A/' A 0\ 

U.0dU(4^ j 


3.219(234) 


11.557(649) 


A OAQ 

U.zUo 


n KAQ / QA^ 

U.oyo^oU J 


2.213(159) 


13.141(758) 


A QC1 

U.ooi 


U.Ooz^zo j 


1.569(131) 


12.354(956) 


A /1AA 


U.00z(^oi j 


1.108(100) 


13.941(1.455) 


A CO /I 

U.DZ4 


A (X 1 1^ / 0A\ 


0.882(74) 


11.980(1.219) 


0.754 


0.484(32) 


0.720(67) 


9.545(1.330) 


1.003 


0.411(42) 


0.460(67) 


8.088(1.810) 


1.122 


0.383(41) 


0.379(52) 


6.238(1.658) 


1.239 


0.368(47) 


0.335(54) 


7.991(2.222) 


1.353 


0.365(53) 


0.312(57) 


7.220(2.114) 


1.465 


0.348(79) 


0.318(83) 


2.968(4.150) 


1.574 


0.315(74) 


0.204(60) 


7.453(3.508) 


1.681 


0.295(89) 


0.196(68) 


4.989(3.480) 


1.888 


0.194(100 


0.112(70) 





TABLE IX: The first column gives the in GeV^, the second 
Og*, the third and the fourth G„ jvA. The errors quoted 
are jackknife errors. 



